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Abstract 

The Galilean transformations of the few-electron atomic wave functions are considered. We 
discuss the few-electron wave functions constructed in the model of independent electrons as well 
as the truly correlated (or highly accurate) wave functions. Results of our analysis are applied to 
determine the probability of formation of the negatively charged tritium/protium ions during the 
nuclear (n, 3 He; t,p)— reaction of the helium-3 atoms with thermal/slow neutrons. 
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The general formulas for the Galilean transformation of the non-relativistic, single-particle 
wave function ty(r,t) are well known since the middle of 1920's (see, e.g., If \I/(r, t) is 
such a wave function written in the reference frames which are at rest, then the corresponding 
wave function $f'(r',t) in the moving frames takes the form 

r 2777 ! 

#'(r-V*,t) = y(r,t)exp[— {-V ■ r + -V 2 t) (1) 

ft 

or 

r 7777, ! 

^'(r, t) = #(r + Vt, *)earj?[— (-V ■ r - - V 2 t)\ (2) 

From Eq.([2]) one easily finds the following transformation formulas for the non-relativistic 
wave function of an arbitrary few-particle system. To avoid analysis of very general quantum 
systems in this study we restrict ourselves to the consideration of few-electron atomic systems 
only. Here and everywhere below by an 'atomic system' we mean the bound few-electron 
system with one very heavy center which also has a positive electric charge Qe. The forces of 
electric attraction between nucleus and electrons bind this system together. The competing 
forces of electric repulsion between electrons decrease the final value of binding energy, but 
they are relatively small. 

The general formulas for the Galilean transformations of the actual few-electron (atomic) 
wave function follow from formulas, Eqs.flT]) - (j2]). For such systems one finds a number of 
advantages to write all formulas in atomic units, where h — l,e — l,m e = 1. In this units 
for the N— electron system we have 

N 

r„ - /V 
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#'(!•!, ...,r N ,t) = *( ri +Vt,...,r N + Vt, t)exp z(-V • - N^-VH)} (3) 



i=i 

This formula is correct for any atomic system which contains iV independent (or quasi- 
independent) electrons. By approximating the actual iV— electron wave functions by the 
trial functions constructed in the model of quasi-independent electrons one can show that 
the formula, Eq.Q, is correct in the general case, i.e. when all electron-electron correlations 
are included. 

In reality, we need the formulas for the Galilean transformations of the iV— electron wave 
functions in the limit £—>■(). This limit corresponds to the sudden approximation 

a, a 

for few-electron atomic systems. For instance, all processes, decays and reactions in atomic 
nuclei proceed significantly faster than usual electron transitions in atoms. Therefore, the 



sudden approximation can be applied to determine the probabilities of the electron-electron 
transitions during nuclear reactions in atoms. In many cases it is important to know the 
probabilities to form various 'final atomic states' after some fast nuclear process in the atomic 
nucleus. The 'incident' atomic state is usualy known. In many cases the newly created nuclei 
(or 'nuclear fragments') are rapidly moving after the nuclear reaction and/or decay. In such 
cases one finds numerous advantages to determine the final state probabilities by using the 
moving frames with the origin located at the atomic nucleus. In these frames all atomic 
electrons suddenly take the speed — V n , where the subscript n means the nucleus. Now, the 
formula for the sudden (Galilean) transformation of the non-relativistic wave function of a 
iV-electron atomic system takes the form 

*'(ri, ...,r N ) = exp(iV n ■ r x + %V n ■ r 2 + . . . + iV n ■ t n )^(y 1 , . . . , r N ) (4) 

This formula is written in the form which can directly be used for an arbitrary iV— electron 
atom with independent and/or quasi-independent electrons. In such cases the wave function 
depends upon the iV electron-nuclear = rj n coordinates (scalars) only. However, all 
actual, few-electron wave functions are truly correlated, i.e. they explicitly depend upon 
both the electron-nuclear r\ = ri n and electron-electron r»j coordinates, which are also 
called the correlation coordinates. Formulas for the Galilean transformations of the electron- 
electron coordinates be obtained from their definitions. Indeed, according to the 
definition of we can write the following identities ry =| r j — r 3 - | = | r$ — r n — (vj — r n ) | = | 
rj — (r n — *V5t) — [rj — (r„ — *V5t)] \ = r^, where 5t is infinitely small. This means that the 
electron-electron coordinates r^- does not depend upon V, i.e. it does not change during the 
sudden motion of the nucleus. It follows from here that the sudden Galilean transformations 
of the truly correlated wave function are also described by the same formula, Eq.(pE]). 

The formula, Eq.fll]), can be applied to determine the probability of formation of the 
tritium 3 H~ ion during the reaction of the 3 He nuclei with slow/thermal neutrons J4] 

3 He + n = 3 H + J H + 0.764 MeV (5) 

in the two-electron helium-3 atom. In the reaction Eq.(jSJ) the notations 3 H and 1 H stand 
for the tritium nucleus (or t nucleus) and protium (or p nucleus). The reaction, Eq.((5]), is 
of great interest for the burning of the high-dense deuterium plasmas [5( (p > 100 g ■ cm~ 3 
(see also pj]). The reaction, Eq.(jS]), in the two-electron 3 He atom and one-electron 3 He + ion 

3 



was considered in our earlier studies [7J1 , ISJ . The cross-section a of this nuclear reaction for 



thermal neutrons with E n m is very large a max ~ 5330 ■ 10~ 24 cm 2 (or 5330 barn) [9]. The 
velocities of the two nuclear fragments formed in the reaction, Eq.([5]), with thermal neutrons 
are v t ~ 1.59632 a.u. and v a ~ 4.78797 a.u. for the tritium and protium nuclei, respectively. 
In this study all particle velocities are given in atomic units, where h = l,m e = l,e = 1 
and the unit of atomic velocity is v e = ac ~ jj= ~ 2.1882661 • 10 8 cm ■ sec^ 1 . Here and 

2 

everywhere below c is the speed of light and a = §- is the dimensionless fine structure 
constant. This 'atomic velocity' v e is the velocity of the Is— electron in the hydrogen atom 
with the infinitely heavy nucleus °°H. It is clear that in atomic units v e — 1. 

Let us evaluate the probabilities of formation of the negatively charged tritium and pro- 
tium ions. In other words, we want to determine the probabilities of formation of the 
two-electron 3 H~ and 1 H~ ions during the nuclear reaction, Eq.([5]), in the two-electron 3 He 
atom. According to the theory of sudden approximations and in respect with Eq.(j3]) such a 
probability of formation of the tritium ion ( 3 H~ or T~) is written in the form Pjj =| A t f | 2 , 
where A/ is the probability amplitude which is written in the form 

Ai f = ($ T -(ri,r 2 ) | #He( r i> r 2)> = ($ T -(ri,r 2 ) | exp(iY t ■ v x + iV t ■ r 2 )# H e(ri, r 2 )> (6) 

where V t is the speed of the tritium nucleus after the reaction, Eq.([5]). In other words, 
the probability amplitude is the overlap integral between the tritium ion and helium-3 wave 
functions, but the wave function of the helium-3 atom must be taken in the moving reference 
frames. 

In the incident and final wave functions we can separate three internal variables r^, 7*31, ri 2 
(or relative coordinates ry =| — | = | ry |) from other six (3 + 3) variables which 
correspond to the translational and rotational degrees of freedom of the whole three-body 
system. Here and everywhere below the notations 1 and 2 mean the electrons, while the 
not at ion/ index 3 stands for the central (heavy) nucleus. In atoms with one heavy nucleus 
the internal coordinates coincide with the interparticle, or relative coordinates. In general, 
the expression for the probability amplitude Aif, Eq.()6]), is reduced to the following form 

Af = / / / 3 ; Z , M( r 3i,r 32 )$ T - (r 32 ,r 31 , r 2 i)exp(iV t ■ r 32 + iV t ■ r 31 ) x 

y lm( t 3i, r 32)^ u e (r 3 2,r 31 ,r 2 i)r 3 2r 31 r2idr 3 2dr 31 dr 2 i (7) 

for the tritium 3 H~ (or T - ) ion. The notations yT M (r^ , r 32 j and y^Zrirw , r 32 ) used in 



this equation designate the corresponding bi-polar harmonics 



lOj (see also They are 



taking care about non-zero angular momenta of the incident and final atomic species. It 
should be mentioned here that the negatively charged hydrogen ion H~ has only one bound 
= 0)— state. Also, in this study we restrict ourselves to the case when the incident 
3 He atom was in its ground l l S(L = 0)— state. In this case all bipolar harmonics in Eq.(J7J) 
equal unity and the probability amplitude, Eq. (J7J), takes the form 

Af = J J J ®T-(r 32 ,r 31 ,r 21 )j (V t ■ r 32 )j (V t • r 31 )^ He (r 3 2 1 r 31 ,r2i)r 32 r 31 r2idr 3 2dr 31 dr2i 
$t- (f*32, r 3i' r 2i)sin(V t ■ r 32 )sin{V t ■ r 3 i)^ He (r 32 , r 31 , r2i)r 12 dr 3 2dr 31 dr 21 (8) 
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V? 

where V t is the speed of the tritium nucleus after the nuclear reaction in the 3 He atom. 

The wave functions of the ground 1 1 S'(L = 0)— states in the two-electron H~ ion and He 
atom are usually approximated with the use of highly accurate variatonal expansion written 
in the relative/perimetric coordinates r 32 , r 3 i and r 2 \ or Ui,u 2 , u 3 (more details can be found, 



e.g., 



12J). The most advanced of such expansions is the exponential variational expansion 
in the relative coordinates. It takes the following form (for the bound S(L = 0)— states in 
the two-electron systems): 

1 ~ N 

ip(r 3 2,r 31 ,r 2 i) = -^[1 + (-1)^12] C l ex'p{-a i r 32 - /3<r 31 - 7^2i)x(l, 2) (9) 
V2 i=x 

where Q are the linear variational coefficients, P\ 2 is the permutation of the two identi- 
cal particles (electrons 1 and 2) and N is the total number of terms in the trial function 
"0( r 32i r 3ii r 2i) which is an accurate approximation of the actual wave function ^f(r 32 , r 3 i,r 21 ). 
In Eq.()9]) the notation %(1, 2) stands for the two-electron spin function. For the singlet states 
X(l,2) = - (5a) and « = in Eq.©. The total energies obtained for the ground 

l 1 S— states of the H~ ion and He atom with the use of Eq.Q can be found in Table I. The 
wave functions, Eq.([9]), are used in calculations of the probability to form the bound tritium 
ion 3 H~. 

For the tritium ion 3 H~ the probability amplitude Aif is written as the double sum of 
the following three-particle integrals 

£$5?i(a, P, 7; V t ) = ^ / / / exp(-ar 32 - f3r 31 - -yr 21 )sin{V t ■ r 32 ) x 

sin(V t ■ r 31 )r 12 dr 32 dr 31 dr 2 i (10) 

Theory of these integrals was developed in [l3|. In particular, it was shown in [13| that such 
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an integral is reduced to the following double sum (here we apply the Cauchy formula) 

OO / 1 \KT/2k K 

BffiCM.TjV) = £ Vtw^CWw(«.A7) (11) 
K =o \ lK + Z J ! ^=0 

where is the binomial coefficient, i.e. the number of combinations from n by k (n > k), 
and Tk-i- n {a, b, c) is the basic three-particle integral defined in Q. This formula allows one 
to determine the probability to form the bound T~ (or 3 H~) and 1 H~ ions during the nuclear 
reaction, Eq.(j5j), in the two-electron 3 He atom. For instance, by using the approximate one- 
term wave functions for the ground state in the helium atom and hydrogen ion given in 
Table II of Ref.j3] we have found that the probability to form the bound T~ (or 3 H _ ) ion in 
the reaction Eq.© is « 0.77048798 % (probability amplitude is « 0.87777445). Such a large 
probability of the T~ ion formation means that these ions formed in the reaction, Eq.flHJ), can 
easily be dected in modern experiments. Analogous probability for the negatively charged 
protium ion is only « 2.391074-10 -5 % (probability amplitude is « 0.48898613-10" 3 ), i.e. it is 



significantly smaller. This illustrates a very sharp dependence o 



the final state probabilities 



upon the final velocity of the atomic particle (for V > 1) Q], {7]. 

It is interesting to note that we can also use the formula, Eq. tfTTT) . in the case V = 0. 
This case corresponds to the f3~ decay of the 3 H~ ion into the two-electron 3 He atom. 
The corresponding probability obtained with our one-term wave functions is ~ 23.893045 
% (probability amplitude is ~ 0.4880502483). These amplitude and final probability are 
very close to our earlier prediction made in 1998. Analogous calculations with the use of 
three-term variational wave functions (with the carefuly optimized non-linear parameters) 
for the H~ ion and He atom gives the following probabilities: 22.065318 % decay), 
0.879736 % ( 3 H~ ion formation) and 5.72186- 10~ 5 % ( 1 H~ ion formation). The corresponding 
variational three-term energies for the °°H _ ion and °°He atom are -0.527665299440825 a.u. 
and -2.90363883557241 a.u., respectively (compare with the energies from Table I). 

Note that our variational few-term wave functions for the both atomic systems (H~) ion 
and He atom) with carefully optimized non-linear parameters are the best functions in their 
class. Nevertheless, it is very interesting to check our predictions by using variational wave 
functions for the H~ ion and He atom with significantly larger number of terms, e.g., the trial 
wave functions from Table I. In actual calculations with such wave functions we have found 
an additional problem related with the use of the non-orthogonal basis sets in calculations of 
the overlap integrals which include the two different wave functions. This problem was never 



discovered in earlier studies where different expectation values were computed for the same 
systems. In such cases the bound state wave functions were exactly identical for the 'incident' 
and 'final' state. The situation changes drastically if we consider the overlap integrals of two 
different wave functions written in the non-orthogonal basis sets. To illustrate the problem, 
let us suppose that we have changed the order of basis vectors in the second wave function in 
the overlap. It does not change the norm of this (i.e. second) wave function, but numerical 
value of the overlap integral can be changed substantially, since the transformation which 
connects these two non-orhtogonal basis sets is not unitary (or orthogonal). Briefly, this 
means that we cannot predict the right order of basis functions (in the two approximate 
wave functions) which is only appropriate for calculations of the final state probabilities. At 
this moment we are trying to solve this intersting problem and develop the new, reliable 
method for accurate computations of the final state probabilities. 
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TABLE I: The total energies E of the ground l l S{L = 0)— states in the negatively charged hydrogen 
ion °°H~ and °°He atom (in atomic units). K is the total number of basis functions used. 

K E(°°H.-) £(°°He) 

3500 -0.527751 016544 377196 590213 -2.903724 377034 119598 030965 

3700 -0.527751 016544 377196 590333 -2.903724 377034 119598 030983 

3840 -0.527751 016544 377196 590389 -2.903724 377034 119598 030995 

4000 -0.527751 016544 377196 590446 -2.903724 377034 119598 031033 



9 



